%           File taxi.m

% X = taxiType = {0 = blu, 1 = verde}
% Y = Colore = {0 = blu, 1 = verde}

clear; 
format compact;

N = 1000000

seed = 0
rng(seed);
ProbTaxiBlu = 0.85 % Probabilità taxi verde
PerformanceTestimone = 0.20 % Probabilità di sbagliare colore
TaxiBlu = 0; TaxiVerde = 0; Corretto = 0; Sbagliato = 0;

for i=1:N
    indice = i
    s1 = rand;      % Taxi blue / green
    if (s1 >= ProbTaxiBlu) 
        taxiType = 1;
        TaxiBlu=TaxiBlu+1;
    else 
        taxiType = 0;
        TaxiVerde=TaxiVerde+1;
    end;
    s2=rand;
    if (s2 >= PerformanceTestimone)
        colore = taxiType;
    else
        if (taxiType == 1) colore = 0; end;
        if (taxiType == 0) colore = 1; end;
    end;
    if (taxiType == 1) && (colore == 1) 
        Corretto=Corretto+1;
    end;
    if (taxiType == 0) && (colore == 1)
        Sbagliato=Sbagliato+1;   % Blu ma era verde
    end;
end;
Attendibilita = Corretto / (Corretto + Sbagliato)

        